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We introduce and study an exactly solvable model of several species of fermions in which 
particles interact pairwise through a mutual magnetic field; the interaction operates only 
between particles belonging to different species. After an unitary transformation, the 
model reduces to one in which each particle sees a magnetic field which depends on the 
total numbers of particles of all the other species; this may be viewed as the mean-field 
model for a class of anyonic theories. Our model is invariant under charge conjugation 
C and the product PT (parity and time reversal). For the special case of two species, 
we examine various properties of this system, such as the Hall conductivity, the wave 
function overlap arising from the transfer of one particle from one species to another, and 
the one-particle off-diagonal density matrix. Our model is a generalization of a recently 
introduced solvable model in one dimension. 



PACS number: 71.10.Pm, 71.27.+a 



Exactly solvable models of interacting particles have often been very useful in illustrating some general 
concepts in many-body physics. While there is a large variety of such models available in one dimension, 
many of which fall into the class of Tomonaga-Luttinger liquidscl, there are few models known in two 
dimensions which are completely solvable. In this paper, we introduce and study a model of several 
species of fermions which interact with each other through a magnetic field term which depends on the 
coordinates of pairs of particles belonging to two different species. The model can be solved by a unitary 
transformation which reduces it to a model of fermions in a magnetic field which depends on the total 
numbers of fermions belonging to the other species. Our model is a direct generalization of the recent 
reinterpretation of the well-known model of Luttinger in one dimensions. The one-dimensional model 
also has pairwise "gauge" interactions depending on the coordinates of the particles; the model is exactly 
solvable because the interactions can be unitarily gauged away at the cost of modifying the boundary 
conditions in a non-trivial way. As we will see, in our two-dimensional model the interactions cannot be 
gauged away in the bulk of the system; the unitary transformation leaves behind a static magnetic field. 

Let us consider v species of fermions in two dimensions (say, the x — y plane), with the charge and 
number of fermions of type a being denoted by q a and N a respectively. The coordinates of the particles 
will be denoted by r*j a , where 1 < i < N a and 1 < a < v. We will consider the Hamiltonian 
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where c is the velocity of light. Bqz is an external magnetic field pointing in a direction perpendicular to 
the two-dimensional plane; we have chosen the symmetric gauge for its vector potential Ai tCt in order to 
explicitly maintain invariance under rotations of the plane. The other vector potential Ai >a arises from 
two-body interactions; it will be taken to have the following form which is natural in two dimcnsionsB, 

A.c* = Va ^2 £,a0 Z X ( T^ a - fj t fi ) , (2) 
30 

where r\ a and are some constants to be fixed below. Note that the Hamiltonian (|l|) is invariant under 
translations in the plane. 

We may now perform an unitary transformation on the Hamiltonian of the form 



h = u nu- 1 



U = exp 



a<{3 i.j 



where q is the charge of an electron. (We note that the phase factor in U only depends on the total 
coordinates R a = r^ Q of the various species of fermions). This gives the transformed Hamiltonian 
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provided that 



iafj — - £,0a , 

and q a ila = Q f° r a ll ■ (5) 



The antisymmetry of implies that the two-particle magnetic interaction can only act between particles 
belonging to two different species. 

It is interesting to consider the effects of some discrete symmetries such as time reversal (T), parity 
(P) and charge conjugation (C). Let us first set the external magnetic field Bq — 0. Under T, the wave 
functions and factors of i are complex conjugated (thus, the momentum operators p^ a — > —Pi,a) and the 
time coordinate t — ► — t; the space coordinates x,y and the various parameters q a ,ri a and remain 
unchanged. Under P, one of the space coordinates, say, x — > —x, while y, t and all the parameters 
remain unchanged. We therefore see that the model is not invariant under P and T separately, but it is 
invariant under the combined operation PT . Under charge conjugation, we demand that q a — > — q a and 
£a/3 - * — £a/3, while r] a and the space-time coordinates remain unchanged; thus the model is invariant 
under C and therefore under CPT. Finally, if the external magnetic field Bq is nonzero, the model is 
again invariant under C and PT, but not under P and T separately; this is because a magnetic field 
(which must be produced by some external currents) changes sign under C, P and T separately. 

It may be useful to point out here that our model has some resemblance to the mean field theory of 
several species of anyons. In the usual theories of anyons, the wave function is assumed to pick up a phase 
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9ij whenever particle i is taken in an anticlockwise loop around particle j, no matter what the size and 
shape of the loop iso. This is often modeled by treating each particle as a point-like composite of charge 
and magnetic flux; when one particles encircles another, the wave function picks up an Aharonov-Bohm 
phase. In understanding the many-body properties of such a system, a fruitful approach has been to 
begin_arith a mean field theory in which the magnetic flux of each anyon is smeared out over the entire 
planeOu. Thus each particle sees a magnetic field proportional to the average density of particles, which 
is similar to our situation. Of course, the analysis of anyons then goes beyond mean field theory to study 
the fluctuations about the average magnetic field, while our simplified model has no fluctuations. It is 
worth remarking that our model has no counterpart for the most popular anyon model which has only 
one species; we need a a minimum of two species. 

To continue, the total magnetic field seen by a particle of type a in our model is given by B a z, where 

B a = Bo + r] a J2 €<*P N P ■ ( 6 ) 

In order to have a well-defined thermodynamic limit N a — > oo, the £ a( 3 must be taken to scale as l/A, 
where A is the area of the system; thus the magnetic field strengths B a in (^|) remain of order 1 as A — > oo 
with the densities p a — N a /A held fixed. We then expect Landau levels to form for each speciesQ. It is 
well-known that each Landau level has a macroscopic degeneracy equal to A\q a B a \/ (2Trhc). The filling 
fraction of fermions of type a is given by 

2-Khc 

■> a = Pa \„ Ft I ■ (?) 

If f a is not equal to an integer for one or more values of a, the ground state of the system is highly 
degenerate. 

For computational purposes, it is convenient to break this degeneracy in one of two ways. We can 
either add a simple harmonic confining potential to the Hamiltonians (Q) and (^) of the form 

n sh = I Yl rL > (8) 
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and take the limit k —> at the end of the calculation, or we can simply impose a hard wall boundary 
condition at some large radius R. Analytically, it is easier to work with the first method since the problem 
of free particles in a combination of an uniform magnetic field and a simple harmonic confinement is 
exactly solvable as we will now discuss. (Let us drop the species label a in the rest of this paragraph 
and in the next). Since the problem has rotational symmetry, the energies and wave functions are 
specified by two quantum numbers, a radial quantum number n = 0, 1, 2, ... and the angular momentum 
/ = 0, ±1,±2, .... If only a magnetic field is present (with, say, the product qB being positive), the 
single-particle states have energies which only depend on the integer n which counts the number of nodes 
in the radial direction; thus 

E Ut i = huj c ( n + i ) , 

, c = ^. (9) 
mc 

In the lowest Landau level (LLL), n = while I can only take non- negative values; all states have the 
energy Eq^ = Tiuj c /2 independent of /. The normalized wave functions in the LLL are given in terms of 
the complex coordinates z = x + iy and z* = x — iy as 



/ qB \(Z+l)/2 7 J - q ts + i 
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where Z = 0, 1, 2, The amplitudes of these wave functions are peaked on circles of various radii centered 

about the origin r = 0; the radii of these "ring" states are given by r; = y/2lhc/(qB). (If qB is negative, 
the LLL wave functions are given by Eq. ( |To| ) with z replaced by z* . Then the angular momentum only 
takes non-positive values). 

If we now add a weak simple harmonic potential mu> 2 r 2 /2 for all the particles, the energies of the ring 
states in the LLL become 



Ei 



y/u% + 4UJ 2 + ( v^c + 4^ 2 - OJc ) 



(11) 



which increase from the origin outwards as \l\ increases from zero. In the many-particle ground state, 
therefore, the fermions fill up the individual ring states from the origin outwards. In the following 
discussion, we will assume this order of filling in the LLL, without explicitly mentioning the simple 
harmonic confinement which justifies it. 

We will now specialize to the case of two species of fermions to illustrate some properties of our model. 
Let us take the masses equal to m for both species, the charges equal to q\ = q2 = q (thus, rji = r\ 2 = 1), 
and the numbers of particles equal to Ni and N 2 for the two species respectively. We will also set 

£12 = -&! = !, (12) 

where 7 is a number of order 1. After the unitary transformation in (||), the two species see uniform 
magnetic fields equal to 

N 2 

Bi = B + 7 -j 

iVi 

and B 2 = Bo - 7 ~j (13) 

respectively. If the number of particles N± = N 2 , the model is invariant under the exchange of the species 
labels 1 <-» 2 and 7 — * —7; this is in addition to the discrete symmetries C and PT discussed in general 
before. 

One of the properties of interest for such a model is the Hall conductivity. In the absence of impurities 
and any other interactions (such as Coulomb repulsion), what is the Hall conductivity of this system if 
the filling fractions fi and f 2 are both integers? It is fairly easy to see that the answer is 

a xy = [ h sign(B 1 ) + h sign(B 2 ) ] ^ . (14) 

This can be derived from the usual formula for the frequency-dependent conductivity 

= £ r WMWM _ (Q|J„|a)<q|Jx|o) 1 (15) 

xy uo ^ I u - E a + E + iff oj + E a -E + ir) J ' [ ' 

where |0) is the ground state of the many-body system, and the sum over \a) runs over all the excited 
states; r\ is an infinitesimal positive number. The current J is given by the second-quantized expression 
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If we now perform the unitary transformation in (||) on both the current and the states, then (|l^) reduces 
to the conventional expression for the current operator of two species of fermions placed in the magnetic 
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fields given by (|l3|). Eq. (|5|) can then be evaluated in the usual wayErEI; in the zero frequency limit, we 
obtain the expression given in (|l4|). The Hall conductivity will remain unchanged if we make our model 
more realistic by including Coulomb repulsion between the particles. 

Another object of interest in this model is the matrix element of the "hopping" operator 

M(r) = c\{r) c 2 (f) (17) 



between the ground state of the system with {N\,N 2 ) particles and all possible states of the system 
with (N\ + 1,N 2 — 1) particles. [The calculation of this overlap is of interest in connection with the 
"orthogonality" catastrophe which is known to occur in Luttinger liquids in one dimension. It may also 
be useful in the context of a two-layer quantum Hall system in which electrons can hop from one layer 
to the other]. Since our original Hamiltonian (Q) is translation invariant, it is sufficient to compute the 
matrix element of M(0) located at the origin. This simplifies the computation for the following reason. 
In a second quantized form, the annihilation operator for any species is given by 

C (T) = ^ VVK^) c n,l , (18) 
it J 



where the sum runs over all one-particle states (n, I) with wave functions tp n ,i, and c nj ; annihilates a 
fermion in the state (n,l). Since only the zero angular momentum states have non- vanishing wave 
functions at the origin, c(0) gets a contribution from only I = but all possible radial quantum numbers 
n. Thus 

c(0) = ^.o(0) c„, , (19) 



where 



|^„,o(0)| 2 = ^f- (20) 

for all n. (This follows from the normalization of the Laguerre polynomials given in Refs. 7 and 8). We 
can now compute the frequency-dependent hopping function 

M(u) = \{a;N 1 + l,N 2 -l\ M(0) |0; A^, 7V 2 )| 2 2n S(hu - E a + E ) , (21) 



where |0; Ni,N 2 ) denotes the ground state of the system with (Ni,N 2 ) particles, while \a; Ni + 1,N 2 — 1) 
denotes all possible states of the system with (Ni + 1, N 2 — 1) particles. 

For simplicity, let us consider the case in which the filling fractions /i and f 2 are both less than 1, and 
Ni = N 2 = N. Then the ground state \0;N,N) is one which the both the type 1 and type 2 particles 
occupy the LLL states with n = and angular momentum I = 0, 1, 2, N — 1. Upon acting on this 
state with the operator M(0) in (p~7|) , we get a state \a; N + 1, N — 1) in which a type 2 particle has been 
removed from the state (0, 0), and a type 1 particle has been added to the state (n, 0), where due 
to the Pauli exclusion principle. Hence the energy difference is 



E a — Eq — —t- 

mc 



(n + | ) - \ \B 2 \ . (22! 



where n = 1,2,3, .... This gives the locations of the 8- functions on the right hand side of (|2l|). We now 
have to find the weights. We can use (|20| ) to show that for species 2, 
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|<a 2 ;iV-l|c 2 (0)|0;AO| 2 = |^ , (23) 

where \d2',N — 1) represents the state in which a particle of type 2 has been removed from the state (0, 0). 
Similarly, for species 1, we have 

+ 4(0) |0;iV)| 2 = ffl, (24) 

where \a%;N+ 1) represents the state in which a particle of type 1 has been added to the state (n, 0) 
where n > 1. To use the results (^3|) and (|24|) for evaluating the matrix elements in ( |2l] ) , we now perform 
the unitary transformation in (0). At this point, we have to worry about two things. Firstly, the N + 1 
particles of type 1 in the states |o) see a slightly different magnetic field than the N particles of type 1 in 
the state |0), since the number of type 2 particles differ by one in the two cases. However, the difference 
in the two magnetic fields is of order l/A which vanishes in the thermodynamic limit. Thus the wave 
functions in the two cases look almost the same since the magnetic field B appearing in (|l0|) differs only 
slightly in the two cases; so when we use Eq. ( p4[) , it does not matter much if we set N 2 equal to N or 
N — 1 to determine the value of B\ given by Eq. (O) . Secondly, we have to worry about the phase factor 
appearing in U which depends on the total coordinates Rn, u — Y]j Tj, a ) recall Eq. (||). At this point, 
another advantage of locating the hopping operator at the origin r = becomes apparent. Namely, we 
see that the phase factors z ■ Rn+i.i x Rn-i.2 appearing in the states \a) cancels with the phase factor 
z ■ Rn,i x Rn,2 appearing in the state |0), since the two states only differ by the addition or removal 
of particles at the origin; this does not change the total coordinate R a of either species. Putting Eqs. 
H) and ( [22| ) together, we see that the hopping function is given by 



n— 1 



(2-Khc) 



We thus get an infinite sequence of ^-functions with equal weight. 

Finally, let us compute the one-particle off-diagonal density matrix for, say, species 1. We assume again 
that N\ — N 2 — N and both the filling fractions f a are less than 1. We have to evaluate 

. N N 

p{r,r)= j Y\d 2 r itl d 2 f j}2 tp*(r, f 2 ,i, fjv.i; r 1>2 , r 2 , 2 , -, r N<2 ) ^(r 1 , r 2) i, fjv,i; f 1)2 , r 2i2 , f N<2 ), 

J i=2 j=l 

(26) 

where we assume that the particles fill up the states I = 0, 1, 2, N in the LLL. We again perform the 
unitary transformation (j^) . The integrand in ( p6| ) then becomes the product of a phase factor 

exp [ ^i.(f-f') x R 2 } (27) 



(where R 2 = ^3,2 and we have used Eq. (12)), four Van der Monde determinants which typically look 
like rifc<;( Zfc - ct — z i.a) an d its complex conjugate for both the species, and the Gaussian factor 

=>i-w<^+ 2 £^>- wx>i,i- (28) 

1=2 j=i 

Since the Van der Monde determinants are invariant under translations, we can immediately integrate 
over the N — 1 independent relative coordinates (i.e., fk,2 — n, 2 ) of the type 2 particles. The total 
coordinate R 2 of species 2 then remains in the form 



G 



When we integrate over R2, we get a Gaussian of the form exp[— d(r — r*) 2 /A], where d is a number of 
order 1. In the thermodynamic limit, we can set this equal to 1 since we can assume that the separation 
|r — is much smaller than the size of the system. We are now left with only the coordinates r^i, with 
i = 2, 3, N, to integrate over. We finally get 

1=0 

where we have assumed that qBi is positive. If we now take the limit TV — > 00, we find that the off-diagonal 
density matrix is the product of a Gaussian times a phase, 



P (f,n ~ §t ex P - ~ :f - - zz») , (31) 



2TTflC 



QBi. , l2 qB 



4hc Ahc 



which is the usual result for a single species of particles in the LLL. 

To summarize, we have introduced and solved a two-dimensional multi-species fermi system with 
mutual interactions of a particular type. The interaction can be converted via an unitary transformation 
into a static magnetic field whose strength depends on the density of particles. The model is quite simple; 
after all, the exact solvability of a Hamiltonian which has a quadratic form should surprise no one. Yet 
the physics of the model is quite interesting. We end up with a strongly non-fcrmi liquid system; further, 
the elements of the orthogonality catastrophe, i.e., a readjusting of all states in response to the addition 
of a single particle, also carry over from the one-dimensional physics of the Luttinger model. Our model 
may thus serve some purpose in understanding the physics of non-fcrmi liquids in higher dimensions. 

We would like to dedicate this paper to the memory of Heinz Schulz. 
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